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Abstract. Using a form of modified dispersion relations derived in the con- 
text of quantum geometry, we investigate limits set by current observations on 
^ potential corrections to Lorentz invariance. We use a phenomological model 

in which there are separate parameters for photons, leptons, and hadrons. 
Constraints on these parameters are derived using thresholds for the processes 
of photon stability, photon absorption, vacuum Cerenkov radiation, pion sta- 
bility, and the GZK cutoff. Although the allowed region in parameter space 
is tightly constrained, non-vanishing corrections to Lorentz symmetry due to 
quantum geometry are consistent with current astrophysical observations. 
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1. Introduction 



The quantum description of gravitation is arguably the largest gap in our under- 
I standing of fundamental physics. In the last decade, a number of lines of research 

fS| I have offered new insights into the theory which will supersede quantum theory and 

I general relativity. For instance, several approaches predict that space is fundamen- 

r"| ' tally discrete. In one approach, eigenvalues of geometric observables have discrete 

Q^l spectra (see e.g. ^ for a review). While it may not be surprising that the quan- 

tization of curved space yields quantized geometry, it is surprising that present 
^ ' day astronomical observations already limit the extent of quantum geometry effects 

These quantum geometry effects arise from an imprint of discrete underlying 
space on propagating modes. Particles with ultra-high energies interact with struc- 
ture on the smallest possible scales resulting in corrections to Lorentz symmetry. 
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^ . Observations on the TeV scale offer an opportunity to test the extent of these 

quantum geometry induced Lorentz symmetry corrections. 

Lorentz invariance may be the most tested symmetry in Nature. Given the 
wealth of evidence in its support, it may seem obtuse to suggest that there may be 
corrections. However, there are reasons to believe that this may not be an exact 
symmetry not the least of which is the fact that any test of Lorentz invariance 
(necessarily at finite energy) leaves an infinite parameter space untested. Due to 
the non-compact nature of the Lorentz group, exact Lorentz symmetry is untestable. 
There are other reasons to suspect exact Lorentz invariance at all energy scales. For 
instance, ultraviolet divergences in quantum field theory point to new physics at 
high energies. Despite these suggestions, only recently has work in quantum gravity 
yielded concrete proposals on how the symmetry might be modified 

In this paper we explore the consequences of modified dispersion relations which 
are motivated by a study of semiclassical states in loop quantum gravity by Al- 
faro, Morales-Tecotl, and Urrutia We study a model in which the usual 
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dispersion relation of special relativity, ~ + m"^ , is modified by leading order 
quantum geometry corrections of the form nipp^. By the semiclassical analysis, 
the parameter k is expected to be of order unity. It determines the extent of the 
corrections to Lorentz symmetry while the Planck length, £p, sets the scale of the 
effects. Due to the sensitivity of process thresholds, the modifications of particle 
dispersion relations are already tested by current astronomical observations. 

What is more spectacular than merely limiting the possible extent of corrections 
to Lorentz symmetry is the proposal p9| ] that this framework may be robust enough 
to elegantly explain three incongruities, or paradoxes, between standard model pre- 
dictions and observational results: (i) Cosmic rays are expected interact with the 
cosmic microwave background (CMB) producing pions and introducing an upper 
limit on the observed energy of particles of cosmological origin. Known as the GZK 
cutoff (2O[0 this upper limit has not been observed. About twenty events at sig- 
nificantly higher energies have been reported |2^, ^ . Modified dispersion relations 
can raise the GZK cutoff |^. (ii) Ultra high energy photons of cosmic origin are 
also expected to interact with infra red background radiation. According to some 
estimates for the background flux, photons of energy 10 TeV or more should not 
be seen due to background induced pair production ||2^. However, higher energy 
events have been reported |p5| , |2^ . Corrections to Lorentz invariance provides one 
explanation for this apparent paradox . (iii) Observations of longitudinal devel- 
opment in extensive air showers of high energy hadronic particles are apparently 
inconsistent with predictions ||2^. As proposed in Ref. |2^, one possible explana- 
tion is that high energy neutral pions become stable. This may also be explained 
using modified dispersion relations . 

We calculate thresholds for processes involving photons, leptons and hadrons 
and, with observational limits, constrain the values for the k parameters thereby 
confining the extent of quantum geometry corrections. In more detail, in the next 
section we summarize the results of Alfaro et. al. ||l6|,|l^. Specifying only general 
properties of a semiclassical state, such as flatness above a characteristic scale L, 
the authors find that, in an analysis of particle propagation, photon and neu- 
trino p^ ] dispersion relations are modified. In Section ^ we give a brief overview 
of threshold calculations before turning to a number of processes including: pho- 
ton stability, photon non-absorption, vacuum Cerenkov radiation for electrons and 
protons, proton non-absorption, and pion stability. We calculate constraints from 
the threshold calculations to investigate whether it is possi ble t hat observed effects 
may be accounted for by Lorentz symmetry corrections. In 3^ we show that asym- 
metric momentum partitioning, first noticed by Liberati, Jacobson, and Mattingly 
|l3| , dramatically affects the constraints on the dispersion relation modifications. 
The results of Section || are summarized in Table [|. Finally in Section ^ we apply 
the constraints together with current observations to limit the extent of potential 
Lorentz symmetry corrections. We summarize the constraints in the final section 
and in Figures || and |[ We find that present day observations tightly constrain - 
but still leave open ~ the possibility of Lorentz symmetry corrections of this form. 

Particularly close to the present work is the paper by Jacobson, Liberati, and 
Mattingly JlBf in which many of these results were summarized. For the most part 
the present work agrees with this paper where the subject overlaps, although this 
work also includes new threshold calculations and constraints for proton vacuum 
Cerenkov radiation, the GZK threshold, and pion stability. 
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2. Modified dispersion relations 

To determine the action of the Maxwell Hamiltonian operator on quantum geom- 
etry Alfaro et. al. specify only very general conditions for the semiclassical state. 
The idea is to work with a class of states which satisfy the following conditions: 
(i) The state is "peaked" on flat and continuous geometry when probed on length 
scales larger than a characteristic scale L, L ^ £p. (ii) On length scales larger than 
the characteristic length the state is "peaked" on the classical Maxwell U{1) con- 
nection, (iii) The expectation values of operators are assumed to be well-defined. 
In addition it is assumed that the geometric corrections to the expectations values 
may be expanded in powers of the ratio of the physical length scales, L and £p. 

States peaked on geometry and the geometric connection are the natural ex- 
pectation for semiclassical or coherent states which model flat space. The work of 
Alfaro et. al. is a forerunner for the detailed analysis of semiclassical states. As the 
specification of semiclassical states becomes more precise, and work is well under 
way [^-^5|, we can expect to check these initial results. 

Expanding the quantum Maxwell Hamiltonian on their states Alfaro et. al. find 
that for massive neutrinos the particle dispersion relations are modified with fl^ 

El=p^ + {k± (3YpP^ (1) 

where p is the magnitude of the 3-momentum, k and [3 are of order unity and 
P is helicity-dependent. The result is computed for a superposition of helicity 
eigenstates. 

For photons the authors find that the quantum geometry effects are given, to 
leading order in £pfc, by 

where T parameterizes the scaling of the semiclassical state expectation value with 
respect to the gravitational connection. As before, the parameters a and (3 are of 
order unity with f3 parameterizing helicity dependent corrections. [| For the purposes 
of exploring the leading order helicity-independent effects in this model of quantum 
geometry, we take T = —1/2 and /3 = 0. Helicity effects have been investigated in 



Ref. 1 36 1 and may be treated using similar methods used in the present work. 

To place all the particles we consider on the same footing we take a phenomeno- 
logical approach in which the modification to the usual dispersion is parameterized 
by a parameter for particle species a. Thus, the effect is simply modeled by the 
modified dispersion relation 

El^pl + ml + Kjppl (3) 

In this equation p is the magnitude of the particle 3-momentum, a quantity greater 
than 0.0 For each particle species, the parameter Ka is, from the semiclassical 
analysis, expected to be of order 1 and can take either sign. Of particular interest 
to determine whether positive allowed for massive particles. In this case, 

massive particles would be able to propagate faster than the low energy speed of 
light, c, which has interesting repercussions for causality. Alternately, for either sign 



^The specification for the class of semiclassical states is more general in ]l7| than in [|l6t . With 
the same specification as the photon case, the neutrino dispersion relation would have the 
same form as (0) jis) . 

Since Eq. M) is in terms of the magnitude, the relation is invariant under time-reversal. The 
modification could also be viewed as a nlpp^E term. 
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of Ka, the particles might be able to propagate faster than high energy photons as 

long as Ka > K^y. 

In terms of bookkeeping, we use £p := y^AnhG/c^ throughout. We set the 
low energy speed of hght to one ("c is c is 1") and fi to one. The scale of the 
model becomes £p ^ 3 x 10~^^ eV^^. The effects of the corrections terms become 
significant when the correction term is of the same magnitude as the mass term, 
i.e. when Pcrit ~ {m^/ipY^^ ^ 10^'^, 10"'^^, and 10^^ eV for electrons, pions, and 
protons, respectively.^ 

The dispersion relation of Eq. (^) is not valid for all momenta. Since the class 
of semiclassical states identified by Alfaro et. al. only reproduces flat geometry on 
scales larger than L, we can only use the dispersion relation within this approxima- 
tion; particles we investigate cannot have wavelengths shorter than L. Thus, the 
momenta are restricted by Pa ^ We remain well within this restriction with 

ipPa < 10-^ 

We call the modification of the dispersion relations "quantum geometry effects" 
since the background quantum geometry only provides the scale on which new non- 
linear terms enter the equations of the motion. In addition the quantum geometry of 
the semiclassical state enjoys no reaction due to propagating modes on the geometry. 
It is clear that beyond the realm of the semiclassical model reaction of the particle 
on the geometry must be taken into account and we should expect new behavior. 

Throughout this discussion we assume that the momenta, energies, and lengths 
are compared in one inertial frame. Indeed, if Lorentz symmetry is inexact then 



there is presumably a preferred frame (however see, for instance, Refs. |39|-41 for 
modifications of Lorentz symmetry which do not select a preferred frame). For the 
purposes of this phcnomenological study, we work in the preferred reference frame 
in which the cosmic microwave background radiation is isotropic. 

Within the context of perturbative quantum field theory it is well known that 
modified dispersion relations can yield causality violations and breakdowns in per- 
turbative quantization (see, for instance, Since we investigate particles for 
which the correction terms play a significant role, these issues are important. How- 
ever, in the following we work in a single inertial frame and do not have the trans- 
formations between inertial frames in this new context. So we leave investigations 
of these issues to future work. It would be particularly interesting to precisely char- 
acterize modified dispersion relations which are consistent within the low energy 
perturbative framework. 

3. Particle Processes 

On account of the modification in the dispersion relation of Eq. (H) quantum 
geometry effects are potentially observable. This is due to the sensitivity of process 
thresholds to the correction term for ultra-high energy particles. We outline the 
framework of the threshold calculations then present the details of each of the 
processes in turn. 

3.1. Threshold Kinematics. Suppose we have the 2 particle interaction a -I- 5 — > 
c + d. We assume that energy is conserved, Ea + Eh = Ec + E^, and 3-momentum 



•^Incidentally, this shows why high precision tests of Lorentz invariance such as in Refs. 
and lis) do not offer strong limits on the ft parameters. The energies are too low. However, if the 
characteristic length scale of the semiclassical state were fixed at much larger scale - say nuclear 
- then these high precision tests would offer limits on the analogous parameters of such a model. 
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is conserved, Pa + P6 = Pc + Pd- The only unusual bit is in the modified disper- 
sion relations. Although the modified threshold calculations are a straightforward 
application of momentum conservation, energy conservation and the modified dis- 
persion relations, the calculations are not without surprises. In fact, even when 
the outgoing particles have the same mass, it may not be energetically favorable to 
partition the momentum symmetrically. This asymmetry was, to our knowledge, 
first observed by Liberati, Jacobson, and Mattingly p^ . 

We almost exclusively use the leading order approximation of Eq. (|^) 

E^p+ — + -Kiy. (4) 

It is clear that this approximation applies for high energy particles only, when 
p ^ m so that m <^ p <^ l/^p- 

The interaction geometry may be strongly affected by the correction term. At 
ultra-high momenta the correction term dominates and, if the sign of k is negative 
then the energy of the outgoing particles may be reduced by including transverse 
momenta. However in the regimes we consider, using 4-momentum conservation 
and these modified dispersion relations it is not hard to see that the threshold 
interaction geometry is what one would expect: the incoming particles' momenta 
are antiparallel and the outgoing particles' momenta are parallel. Of course, the 
momenta must be sub-Planckian. We also emphasize that, in these dispersion rela- 
tions, the K-parameters are of order one. Since the framework of the semiclassical 
states requires that the momenta be far below the Planck scale, we employ this 
interaction geometry throughout the remainder of this paper. 

We use the modified dispersion relation Eq. (jj) in the conservation of energy. 
With energy and momentum conservation, we find the kinematic constraint 

7^ + T^^JpPI + 7^ + T^^^b^pPb = 77^ + T^^c^pPI + 77^ + T^^dpPl (5) 
Ipa I Ipb ^ ^Pc ^ ^Pd ^ 

By momentum conservation Pd ^ Pi ^ Pc with pi '■= Pa + Pb being the available 
incoming momentum. Thus the right hand side of Eq. may be expressed as a 
function of Pc only. 

To locate the threshold for one of the incoming particles (typically, the other 
momenta is known) , we wish to find the minimum final energy as a function of the 
momentum of the outgoing particle c. This may be accomplished by differentiating 
Eq. (||) and finding the roots. In the general two-channel case, this is given by 

9 9 
- -TT^ + Kjpp, + ^ + Kdp{p^ - Pc) (6) 

2p2 2{p,-pcy 

which is a quintic in pc- With the physical root(s) in hand one can return to Eq. 

to find the threshold momentum. However, in the individual processes there 
are better methods than attempting to directly solve this quintic. Perhaps the best 
way to see this is in the case of photon stability. As this case is straightforward 
and yet gives asymmetric partitioning we present it in its entirety. 

3.2. Photon Stability: 7 7^ e+ -t-e~. In special relativity photon decay is forbid- 
den - a simple consequence of the energy and momentum conservation. However, 
taking quantum geometry effects into account it is energetically favorable for ultra- 
high energy photons to decay. Recent observations of multi-TeV photons [^,|2^, ^ 
provide potential threshold values for the quantum geometry induced decay. 
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Applying the threshold framework above, we immediately see that there are 
considerable simplifications in Eq. (^). In fact, the energy- momentum constraint 
becomes 



2' 1 1 

5^ = "I 

\Pe+ Pe 



+ nJp{pl++Pl-)- (7) 



Before differentiating, notice that for negative K-y, the electron K-parameter must 
also be negative. This opens up the possibility that a simple partition of 

incoming momenta may not be energetically favorable. 

This is already clear from the modified dispersion relations for negative Ka- As 
the momentum of the particle increases the energy increases. However, when the 
momentum is near Pcrit the energy increases less rapidly]^ This curve strongly 
affects the nature of threshold calculations: At high momentum the slope of the 
energy steadily decreases so it becomes energetically favorable to partition the mo- 
mentum of the outgoing particles asymmetrically. By asymmetrically partitioning 
momentum one can only slightly increase the energy of the high momentum particle 
while drastically reducing the energy of the lower momentum particle leading to a 
lower total energy (for more extensive discussion see Ref. Q). 

With this observation we re-express the outgoing momenta as p^- = Po — ^ 
and Pe+ = Po + A with po being half the available momentum, po '■— P7/2. The 
asymmetry factor A can be no larger than p^; —p^ < A < po- The relation Eq. (|^) 
becomes 

2'^7^pPo = + ^e^M + A^). (8) 

To find the threshold energy, or minimum energy, for this process we minimize this 
equation. This yields a simple quartic for A which has two physical roots: A* = 
for the symmetric and 



for the asymmetric partitioning, respectively. The second root is only real for 
negative Ate- Since the momenta asymmetry is bounded by po there is an upper 
bound on Kg for this root: 

In addition, when the threshold energy is found using asymmetric root another con- 
dition is required to ensure that all the momenta are non-negative. This condition 

To find the regions of K-parameter space where the two roots impose constraints, 
it is worth returning to Eq. (^. As observed above, if is negative then Ke 
is as well. It is then energetically advantageous to asymmetrically partition the 
momentum, i.e. to use the non-zero root A^. Thus, when both < and the 



*At high momentum, at the limit of the semiclassical approximation, the energy reaches a max- 
imum. Understanding what occurs at these energies requires detailed knowledge of the dynamics 
in addition to the semiclassical state of quantum gravity. 
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inequality of ( [lO| ) is satisfied the asymmetric configuration is applicable. With the 
two roots and the energy- momentum constraint, the momenta at threshold are 

-,1/3 



P7- 
and 



8mi 



for K-y > 



(11) 



1/3 



for 



< < 



where the two thresholds are valid in the regions identified by the inequalities. 

Given a threshold value for the process, these equations give relationships 
between Kg and Kj . Since the photon is stable if the incoming energy is less than 
the threshold energy, we can assign inequalities to the constraints on He and k^. 
Thus, the photon is stable at the threshold p-y^ if 

4™-^ 



1 

2 



f-ppl 



< for > and 



(12) 



< for Ke < K-y < 0. 



The regions of applicability of the two constraints are inherited from the momenta 
thresholds. These constraints, together with the regions of applicability, determine 
allowed regions in the k parameter space. 

The highest energy observed gamma ray, 50 TeV, originated from the relatively 
local Crab nebula . Using this energy for p^^ we plot the constraints of (|l^) in 
Fig. (|^). The constraint inequalities ensure that the photon does not decay below 
this threshold. The shading represents regions in — space which are ruled out 
by the constraints of (|l2|). The smooth transition from the linear to the quadratic 
constraints on the k parameters occurs when changes sign. The resulting curve 
forms the upper bound on the allowed region. 

There is another process which has a closely related calculation. This is the 
apparent paradoxical observations of TeV photons from distant sources. 

3.3. Photon non-absorption: 7 + 7//? 7^ e+ + e^. Due to the presence of back- 
ground radiation, high-energy photons are expected to be absorbed in pair creation 
by the far infra-red (IR) background. For photons originating from distances greater 
than ^ 100 Mpc, the 0.025 eV background provides a cutoff in the spectrum at 
^ 10 TeV. However, recent observations of multi-TeV photons have challenged this 
expectation [^^. Located at redshift ^ 0.03 the active galactic nuclei Markarian 
421 and Markarian 501 have produced flares with multi-TeV photons. Observations 
have included 17 TeV |2^ and 24 Tev ||2^ events. It has been suggested that the 
apparent overabundance of high energy photons may be due to to corrections to 
Lorentz invariance along the lines of Eq. (|) 0. 

The process kinematics is identical with the case of photon stability of Eq. (||) 
except the presence of the IR photon. Letting the energy of this low-energy photon 
be e, conservation of momentum for this process becomes p^ — e — Pe+ + Pe- and 
equation for photon absorption becomes identical to Eq. (^) with the substitution 
ll 



4e 
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Figure 1 . Using the thresholds for photon stability and a thresh- 
old value of 50 TeV, the two constraints are plotted. The two curves 
smoothly intersect at the transition from symmetric partitioning, 
the diagonal line, to asymmetric partitioning, the quadratic func- 
tion in the third quadrant. The vertical dashed line is the condition 
of Eq. (|lO|). This identifies the upper branch of the asymmetric 
constraint as the physical solution. The region ruled out by photon 
stability up to 50 TeV is indicated with shading. 



Thus, the calculation of the threshold is the same as for photon stability and result- 
ing K inequalities are identical with the substitution ^ k'^. Because of the factor 
of Po in the definition of kI^, it is not as straightforward to write simple expressions 
for p^, . Our main interest is on the allowed regions in n parameter space so we are 
satisfied with the equations which relate Kg and Kj for the allowed region. These 
are presented in Table ||. Plots of the allowed region are similar in form to the ones 
for photon stability. 

3.4. Pion Stability: tx ^ Within the context of ordinary particle kinemat- 

ics, the neutral pion has a lifetime of 8 x 10^^^ s. As suggested by Amelino-Camelia 



19 1, the modified dispersion relations can have the effect of making ultra- high en- 
ergy neutral pions stable. This may explain some inconsistencies between observed 
and expected patterns in cosmic ray showers. 

Normally the pion decays because the rest mass can be converted into energy 
in the form of photons. With the modified dispersion relations, it is possible that, 
above a threshold energy, it is no longer kinematically possible for the pion to decay. 

The threshold calculation is similar to photon stability. Proceeding as before we 
asymmetrically partition the available momentum for the two photons, Po — A and 
Po + A. The resulting energy- momentum constraint is 

2 

5^ + ^K^l^vl = 2At^^p(p2 + A2). (13) 

As is clear from this equation the analysis splits into two cases depending on the 
sign of K-y. For K-y > the minimization requires symmetric partitioning: A* = 0. 
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Solving Eq. ( p^ for the threshold momentum we have 

-,1/3 



Ptt, 



2mi 



for > and > 2k^. 



(14) 



The last inequality ensures that the momentum is positive. 

In the second case, < 0, the symmetric partitioning is in fact the worst choice 
we could make. Equal photon momenta is the maximum of the outgoing energy. 
The local minimum for the energy of the outgoing photons is at the maximum 
possible asymmetry, A* = Po = Pir/2, corresponding to the threshold emission of a 
zero momentum photon. From Eq. (13|), the threshold momentum for this process 



is bounded from below 



> 



1/3 



for < and > 



(15) 



The first inequality for ensures that the energy remains positive while the last 
ensures that the threshold momentum is positive. 

In an identical manner to the previous calculations, Eq. (14]) and ( p^ ) yield 
inequalities among and n^. These are summarized in Table ip 

3.5. Vacuum Cerenkov Radiation: + Cerenkov radiation is observed 

when a charged particle enters a medium with a speed greater than the phase 
velocity of electromagnetic waves in the medium. Vacuum Cerenkov radiation in 
which empty fiat space constitutes the medium is normally forbidden by ordinary 
threshold kinematics. However, Coleman and Glashow suggested that within the 
context of a variable-speed-of-light theory, vacuum Cerenkov radiation could be 
induced when a charged particle's speed exceeds the local speed of light (see 
also [|ll|). Vacuum Cerenkov radiation (VCR) may also be kinematically allowed 
by the modified dispersion relations. This is clear even from a glance at a charged 
particle's speed 



dpa 



1 9 ^ l^a^pPa 

^pI 



(16) 



where a is the charged particle. 

Since both charged particles and photons have the same functional form of the 
dispersion relations, there are four possible cases to investigate. These depend 
on the signs of and Ha' (i) For k,^ > and Kq < no radiation occurs, (ii) 
For K-y > and > vacuum Cerenkov radiation may occur when the charged 
particle's speed exceeds 1, the low energy speed of light. We call this Type I 
VCR. (iii) For < Ka vacuum Cerenkov radiation occurs either when the charged 
particle's speed exceeds the low energy speed of light or when it exceeds the speed 
of the ultra-high energy photon. We call this last case Type II VCR. These cases 
show that VCR could occur in the first, fourth, and part of the third quadrants. 

We present the threshold calculation first a charged particle a. In Section |^ 
we use these results for electrons and protons. The threshold calculations can be 
found using the methods use employed above or, equivalently, using the Cerenkov 
condition Va > v-y. Using the later method we find that for high energy charged 
particles the condition for Type I VCR is 

™2 ll/3 



Pa, 



for Ka > 0. 



(17) 
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This condition corresponds to the momentum at which the particle's speed exceeds 
1, the low energy speed of light. This process the photon emission begins with long 
wavelength photons. For this Type I radiation the spectrum extends from zero 
energy up to = Pa, ■ 

When < Ka, however, the emission process starts with finite energy photons. 
For high energy particles and photons - what we call Type II VCR - the Cerenkov 
condition yields 



Pa, = 



-,1/3 



for Ka > 



3.2 



(18) 

The observational 



This may also be derived using the methods of Section 
constraints may be derived from this threshold. 

When using the velocity relation, one must also ensure that energy is conserved. 
We study the case in which both outgoing particles have high momenta. In the 
leading order approximation of Eq. (^) , the energy-momentum constraint takes on 
the familiar form of the threshold calculations: 



Pa 



+ KaipiPa ~P-/f + ' 



(Pa-P-r) 

After some simplification, this gives a quadratic for 



4- 



(3Ka + Hj)pa 



-P-r 



0. 



(19) 



The requirements that the physical root of this equation, , be real and positive 
produce the threshold value 



Pa, 



[Ha 



1/3 



for Kj < 



-3Ke < 0. 



(20) 



The last inequality arises from the positivity of the maximum energy outgoing 
photon 



Pi, 



(3K(j + Kj) 
2(^a ~t- /^^) 



Pa,- 



(21) 



The two types of VCR have distinct spectra. For Type I, the emission of photons 
starts with very low energy photons and extends up to the particle's energy Pj^^^^ — 
Pa ■ For Type II VCR, the spectrum forms a band which begins at the finite energy 
corresponding to the condition Va{p) — c-y (p) and extends up to the particle's energy 
Pa- The half angles of the Cerenkov cone apertures also differ for the two types of 
radiation. The angle may be simply computed from Eq. (p^). For Type II, 



Oc 



— (1 KjipPy^ '^f^a 



1/2 



(22) 



The angle for Type I is found by setting = in the above equation. 

Using observations of high energy electrons and protons, we can establish limits 
on the the possible values of Kg and Kp. This is done in the final section. 
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3.6. Proton Non- Absorption: p + 7 7^ p + tt'^. Shortly after the cosmic mi- 
crowave background was discovered, Greisen |Q and Zatsepin and Kuzmin 
predicted that the cosmic ray spectrum should have a cutoff. They found that 
cosmic rays constituents including protons and some heavier nuclei arriving from 
cosmological distances would lose energy through interactions with the CMB pho- 
tons. The dominant process, photopion production, provides an effective energy 
cutoff for high energy protons from cosmological sources. In fact the spectrum was 
believed to be bounded from above by 5 x 10^^ eV. However, the observed spec- 
trum extends well above this cutoff to 3 x 10^° eV. While the number of events with 
energies above the GZK cutoff is not large (~ 20), the locations on the sky are not 
inconsistent with the hypothesis that the sources are isotropic. If these events are 
indeed particles originating from cosmological distances then the GZK cutoff must 
be modified. As observed in Q, the modified dispersion relations considered here 
can raise the GZK cutoff. 

The method for computing the threshold energy for photopion production is 
similar to the above process but the details make the problem more complex. The 
different masses of the pion and the proton mean that much of the simplifications 
used above do not work. This returns us to the general two-channel case and the 
resulting quintic for the momentum pc- In addition, when Kp is negative, large 
asymmetries are favored. In fact, the threshold process produces a low momentum 
pion which means that we must use the dispersion relations of Eq. (Q) instead of the 
leading order approximation. Hence the constraint on Kp was found numerically. 

Using threshold values in the range 1 x 10^^ - 8 x 10^° eV we performed a search 
for the value of Kp which separated the parameter space into a region in which the 
high energy proton would be stable and a region where it would be unstable. The 
selection criteria were simply that Ep+e > E'p+Ej^, where the energies are expressed 
using the modified dispersion relations of Eq. (^), and e is the energy of the cosmic 
microwave background photon, 7 x 10^'* eV. Conservation of 3-momentum was 
also required. Since the k parameters for both pion and proton are identical, the 
procedure produced a constraint on Kp as a function of proton threshold energy. 

This function has the familiar shape of the effective potential of a particle in 
a central potential. At energies below the GZK cutoff Kp is positive, effectively 
lowering the threshold. At the GZK cutoff Kp passes through zero. However, 
the dependence of Kp on the threshold energy is not a monotonically decreasing 
function. To see this note that since the magnitude of the proton correction term is 
larger than the magnitude of the correction terms of the products, at energies above 
the GZK cutoff Kp must be negative. At significantly higher energies Kp approaches 
zero from below (as may also be seen with energy-momentum conservation in the 
leading order approximation). The resulting function has a minimum of —7.9 x 
10^^^ at 1.6 X 10^° eV. Thus, measurements of higher energy cosmic rays will 
not restrict the parameter further. These results are discussed further in the next 
section. 

4. Limits on k-parameters 

With all the above processes, one can derive limits on the modifications to the 
dispersion relations. The threshold constraints of Section ^ determine the allowed 
regions in the and parameter spaces. These constraints are summa- 

rized in Table 1^. With the threshold energies in hand one can plot the constraints 
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Process Type 


Constraint 


Applicability 


Photon Stability (7 7^ e + e+) 


Symmetric momenta 
Asymmetric momenta 




K-y > 

^ ^ 


Photon non-absorption (7 + 7/^^ 7^ e + e+) 


Symmetric momenta 
Asymmetric momenta 


(k-v + 7^ Ke)' + < 


ft- > "le 

^'^ ^ (pP-y, 


Pion Stability (tt 7 + 7) 


Symmetric momenta 
Asymmetric momenta 


^ 1 ™^ 


K-y > 2K-n and K-y > 
> and < 


Proton Stability {p + e p + n) 


Numerical Result 


Kp < -8 X 10-16 




Vacuum C'erenkov Radiation {a -/^ a + "f) 


Type 1: Zero-energy photon 
Type 11: Finite-energy photon 
Energy Conservation 


* ™" 

(kq k^) -|- ^^p3 -|- k-j,) < 


K-y > or Kg > 

K-y < 



Table 1 . Summary ol constraints and their regions ol applicability. 



and determine the allowed regions. But before discussing these plots, we summarize 
the observational data which provid e th reshold values for each of the processes. 

We already discussed in Section 3.2 the observation which we use for photon 
stability (the 50 TeV event from the Crab nebula). We now turn to photon non- 
absorption by the far infra red background. 

The unexpected transparency of the universe for high energy photons may be due 
Lorentz symmetry corrections (see, however, Ref. [^). Although direct measure- 
ments from two instruments onboard COBE and detailed calculations have been 
completed, the IR flux is not well-determined |^]. This is largely due to, on the 
observational side, the flux of IR radiation produced by interplanetary dust and, on 
the calculational side, uncertainties in the galactic-evolution models. Given the un- 
certainties in the IR flux, we only plot an example case. The photon non-absorption 
constraint has a complex dependence on energy. Nonetheless, in the relevant region 
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the threshold energy which aUows for the maximum increase in transparency due to 
suppressed pair creation induced by background 0.025 eV IR photons is 15 TeV. In 
the summary plots we use this energy as a representative case. If the unexpected 
transparency is due to Lorentz symmetry corrections then this constraint is the 
lower boundary on the allowed region. 

At shorter wavelengths the spectrum is much better established. In fact, with 
no modifications the universe becomes highly opaque to photons above 100 TeV 
from the interactions with the 2.7K cosmic background radiation |4^. Observations 
of photons above this "hard limit" would necessitate more serious consideration of 
photon non-absorption due to modifications to Lorentz symmetry. 

As shown in Section p.4| modified dispersion relations can make the pion stable 
(see also |Q,|l^]). Recently Antonov et. al. reported on a test of Lorentz invariance 
using longitudinal development of extensive cosmic air showers |2^ ]. They found 
that simulated air showers produced a better fit to the observed depth of the highest 
energy air shower if the pion became stable at high energies. This stability could 
be achieved if pions are stable above pjr. ~ 10^® eV. We use this threshold value in 
the plot. This becomes a very restrictive limit. 

Vacuum Cerenkov radiation provides a strong restriction "from the bottom" in 
Ka — parameter space. We consider limits on the allowed region in the pa- 
rameter space for both electrons and protons. The limit for electrons comes from 
observations of a supernova remnant. The best fit for the X-ray spectrum from the 
supernova remnant SN 1006 arises from assuming that the radiation is synchrotron 
radiation from electrons with energies up to ^ 100 TeV The highest energy 

primary cosmic rays are believed to be protons [ pTf . As we see next the cosmic ray 
spectrum reaches 3 x 10^° eV. We use this limit for VCR for protons. 

Finally, cosmic rays above the GZK cutoff are observed. This cutoff may be 
raised by modified dispersion relations. Using the highest observed cosmic ray, 3 x 
10^" eV for the threshold, we would obtain the constraint Kp < —5.6 x 10~^^. There 
are two interesting features to this result. First, there is minimum at 1.6 x 10^° eV 
in the Kp-energy curve discussed in Section 3.6. To prevent a gap in the cosmic ray 
spectrum between 1.1 x 10^" and 3 x 10^° eV we use the value of Kp at the minimum, 
—8 x 10^^^. Second, even when one accounts for the uncertainties reported in Ref. 
|22| the constraint is not consistent with Kp = 0. Attributing the raising of the GZK 
threshold modifications of the dispersion relations restricts the hadronic parameter 
to negative values and effectively rules out exact Lorentz invariance. 

In Figures (||) and (|^) we identify the allowed regions of and 
rameter space. We include all the above processes, including the more speculative 
processes of the photon non-absorption, proton non-absorption, and pion stability. 
This allows us to check whether a consistent set of parameters exist which simulta- 
neous explain the three paradoxes mentioned in the introduction and are consistent 
with current observations. 

In the electron-photon plot of Figure ^ we include constraints from photon sta- 
bility, vacuum Cerenkov radiation, and photon non-absorption. The allowed region 
is a tightly constrained area mostly lying in the third quadrant and containing the 
origin. Thus, based on these processes it appears that quantum geometry effects 
reduce the energy (and speed) of particles at high momentum. As Fig. ^o) reveals, 
however, there is a region in which positive values of the parameters are allowed. 
The allowed region contains the origin and a tiny sliver of positive Ke corresponding 
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(a.) 




-0.04 



(b.) 



0.04 



Figure 2. Plots of the allowed regions in K,y — Ke parameter 
space. The unshaded portion shows the region in which corrections 
to Lorentz invariance are consistent with current observational lim- 
its. In (a) limits from photon stability, vacuum Cerenkov radia- 
tion, and photon non- absorption are shown. The vertical shading 
corresponds to the photon stability constraint while the horizon- 
tal shading corresponds to vacuum Cerenkov radiation. The con- 
straint labeled / arises from assuming that high energy photons 
are not absorbed in pair creation. If so, this provides a lower limit 
on the allowed region. In (b) the same constraints are shown in a 
region around the origin. The shading is the same as in (a). 



to superluminal electrons (ke < 4.5 x 10^^). The parameter is not as tightly 
confined around the origin and is less than ^ 0.03. 

In the hadron-photon plot of Figure || the limits are obtained from the processes 
pion stability, proton non-absorption, and vacuum Cerenkov radiation. The qua- 
dratic VCR constraints do not show up on this scale due to the high upper limit on 
proton energy, 3 x 10^'^ eV. In this plot the allowed region is a exceedingly narrow 
band again mostly in the third quadrant; the scale is on the order 10"^". There is 
a small region of positive values of the parameters allowed. The band of allowed 
values is constrained by the pion stability from above and VCR from below. 

If these two constraints are borne out by further observations, different modifi- 
cations of Lorentz invariance of this form for hadronic and massless particles would 
be effectively ruled out. It is also interesting that the case Kp = Kj lies nearly on 
the boundary of the allowed region. The difference of this relation and the VCR 
constraint is ^ 10~^^. 

We also note that the allowed region does not contain Kp = 0, the Lorentz 
symmetric case. This is due to the raising of the GZK threshold. It is likely 
that the GZK threshold may be explained by other means. Nonetheless, if the 
GZK cutoff is raised by corrections of this form, then current observations rule 
out exact Lorentz symmetry. So if we take the three conditions of pion stability, 
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Figure 3. The plot of the allowed region in the 
Constraints arising from pion stability, proton non-absorption, and 
vacuum Cerenkov radiation are shown. The vertical shading is due 
to the pion stability constraint while the horizontal shading comes 
from VCR. The constraints due to proton non-absorption and VCR 
appear to be a single vertical line. 



proton non-absorption, and VCR seriously, then the Lorentz symmetry is broken 
and the parameters are extremely finely tuned to each other. From the plot and 
the constraints on k from the semiclassical analysis, we conclude that Kp = 0. 

Finally, it is also clear that nearly all the positive k space is ruled out. In 
particular, the results for massive particles effectively rule out "faster than light" 
behavior. One reason why this is significant is that satisfies the classical notion of 
causality thereby satisfying one of the requirements for consistency in perturbative 
quantum field theory. 

It is important to keep in mind that these are preliminary results. Our strategy 
is to combine all processes to see whether the results are consistent. But it is clear 
that these constraints are not equally precise. For instance the vacuum Cerenkov 
radiation constraints are tentative. The limits for electrons are derived from indirect 
measurements of synchrotron radiation around a single supernova remnant. We 
need more data on high energy electrons and cosmic rays before these processes 
can truly rule out regions of the parameter space. 

Our VCR limit for hadrons could also be called "soft." We assume that the 
highest energy cosmic ray event was due to a proton. As new cosmic ray data 
becomes available we will be able to better judge the validity of this assumption. 
Once the Pierre Auger project ]49| ] is running we should have a much better idea 
of the value of the threshold. Aside from data from more ultra-high energy events, 
this detector should also determine the mass composition of the primary cosmic 
rays. This is absolutely critical for the study of Lorentz symmetry modifications. 

Much better estimates of the far infra red background are needed before we can 
establish the lower limit of the allowed region in Ke — K-y parameter space. Without 
this bound we can not confine the parameters to a finite area. In Ref. [ISl the 
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uncertainty is the background is handled in a different way. Assuming that the 
actual threshold lies between 10 TeV and 20 TeV they confine the allowed area to 
a finite size. 

The pion stability constraint will also be improved with new cosmic ray data. 
The results of Ref . [pTt compare simulated showers with the single, highest energy 
cosmic ray event ]22|]. It will be very interesting to see whether the onset of pion 
stability is observed in other ultra-high energy showers. 

The only constraint which could be considered a hard limit comes from photon 
stability. This sharply restricts the possible superluminal behavior of ultra-high 
energy photons. As higher energy gammas are observed from local sources, this 
constraint will "push down" on the allowed region. 

5. Conclusions 

Beginning with a form of modified dispersion relations, we use exact energy- 
momentum conservation, to derive thresholds for particle processes. The dispersion 
relations we consider are motivated by a class of semiclassical states of loop quantum 
gravity which are classical and flat above a characteristic length scale ij] . We 
work in a model in which there are separate parameters for modifications of the 
dispersion relations for photons, leptons, and hadrons. The process thresholds 
provide a sensitive test of the quantum geometry effects. In some cases, such 
a photon decay and vacuum Cerenkov radiation, normally forbidden processes are 
activated at high energy. In other cases, such a pion decay, processes are suppressed. 

Current TeV-scale observations offer severe restrictions on the nature of these 
effects. There exist allowed regions in both parameter spaces consistent with all 
of the above processes. The results are summarized in Figures || and || and in the 
last section. In particular, including all of the processes the analysis shows that 
the hadronic parameter is effectively equal to the photon parameter and both of 
these are non-zero. Further the lepton parameter is confined to a band in the third 
quadrant as shown in Fig. ||. 
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